Complex  Spatial/Temporal  CFAR 

Z.  Ebrahimian 

Communication  Sciences  Institute,  Department  of  Electrical  Engineering-Systems 
University  of  Southern  California,  Los  Angeles,  California,  90089-2565 
zebrahim@usc.edu 


Abstract 

The  conventional  cell  averaging  Constant  False  Alarm  Rate  (CFAR)  criterion  and  its  variations  work  well  only  in 
strictly  spatially  stationary  environments.  In  non-homogeneous  environments,  clutter  map  (scan-by-scan)  processing  is 
deployed.  The  performance  of  this  method  degrades  in  the  presence  of  slow  targets.  In  this  paper,  a  hybrid  procedure 
for  CFAR  is  proposed,  which  combines  the  advantages  of  both  spatial  and  time  averaging.  The  detection  probability  is 
derived  and  the  related  plots  are  given  for  different  values  of  L ,  the  number  of  persistence  scans.  A  method  is  presented 
to  choose  the  parameter  of  the  hybrid  CFAR  to  have  the  lowest  self-masking  effect. 


I.  Introduction 

The  returned  signal  in  a  radar  receiver  is  composed  of  target  reflections,  signal  reflection  from  clutter 
points,  and  thermal  noise  in  the  receiver  input.  With  a  constant  threshold  in  a  clutter  environment, 
false  alarm  probability  is  extremely  sensitive  to  variations  in  average  power  of  the  clutter  signal. 
Although  it  is  possible  to  decrease  the  false  alarm  rate  by  choosing  a  high  threshold,  this  will  decrease 
system  sensitivity  to  the  target  in  low  clutter  regions,  and  reduce  detection  probability.  Accordingly 
a  constant  false  alarm  rate  (CFAR)  approach  is  used  to  adaptively  select  the  detection  threshold.  In 
these  types  of  detectors,  in  order  to  declare  a  target,  one  checks  for  presence  of  the  target  signal,  in  the 
cell  under  test,  which  is  strong  enough  for  the  reference  cells.  Three  important  ways  for  determining 
the  threshold  in  an  adaptive  manner  are  parametric  CFAR,  non-parametric  CFAR,  and  Clutter-Map. 
In  the  two  first  cases,  it  is  assumed  that  the  samples  in  the  adjacent  cells  to  the  test  cell  (which  are 
known  as  reference  cells  and  can  be  expanded  in  the  range,  azimuth,  or  doppler  domain)  are  iid,  and 
a  function  of  those  will  be  used  as  threshold.  The  main  assumption  in  parametric  CFAR  is  that  the 
probability  density  of  interference  is  known  but  its  parameters  are  not  identified.  Reference  window 
cells  are  then  used  to  estimate  the  unknown  parameters  and  the  threshold  is  determined  based  on  these. 
Fig.  1  shows  a  simple  CA-CFAR  (Cell  Averaging)  of  this  group.  This  method  is  predicated  upon  the 
presence  of  clutter  with  stationary  statistics.  As  such,  in  the  case  of  moving  rainstorms,  or  jamming 
and  interference  with  signals  of  other  radars,  the  false  alarm  rate  may  rapidly  increase.  Also  these 
methods  show  a  considerable  drop  in  quality  in  the  clutter  edge  or  in  the  case  of  multiple  interfering 
targets.  Clutter-MAP  was  produced  to  counter  with  non-homogenous  clutter.  In  these  methods,  the 
signal  level  in  the  cell  under  test,  over  successive  scans,  is  used  to  estimate  the  interference  level  instead 
of  averaging  on  the  adjacent  cells.  This  method  is  not  affected  by  non-homogenous  backgrounds  and 
increases  sub-clutter  visibility  ([1],  [2]).  In  non-homogenous  environments,  CM-CFAR  (Clutter  Map- 
CFAR)  shows  much  better  performance  than  CA-CFAR.  In  CM-CFAR,  the  assumption  is  that  the 
clutter  samples  are  independent  over  scan  by  scan. 


Fig.  1.  CA-CFAR  processor 


Fig.  2.  CM/L-CFAR  block  diagram 


In  spite  of  the  above  mentioned  advantages  for  methods  based  on  clutter-map,  their  capabilities  are 
limited  by  high  amounts  of  memory  requirement,  transition  response,  and  velocity  response.  Detection 
loss  for  slow  targets  is  increased  by  velocity  response  limitation  and  thus  causes  the  self-masking 
problem.  To  improve  these  drawbacks,  complex  methods  of  spatial  and  time-averaging  CFAR  are  used. 
One  of  these  complex  procedures  was  proposed  in  [3]  (Fig.  2).  In  this  method,  sample  cancellation  by 
the  Ll-Filter  [4]  before  the  clutter  map  is  a  very  useful  way  to  decrease  self-masking  and  interfering 
target  masking. 

In  this  paper,  we  propose  a  new  complex  method  of  spatial  and  temporal  CFAR  (see  Fig.  4).  In 
this  method  some  percentage  of  the  detection  threshold  is  obtained  from  CM-CFAR  and  the  rest 
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from  temporal  CFAR  (thus  a  +  f3  =  1).  This  preserves  the  CM-CFAR  immunity  against  spatial 
non-homogeneities  and  simultaneously  exploits  the  good  velocity  response  of  spatial  CFAR. 

The  rest  of  the  paper  is  organized  as  follows.  First  the  theoretical  analysis  of  the  new  method  is  given. 
Then,  the  self-masking  effect  on  the  system  performance  is  investigated.  The  experimental  results  is 
described  at  the  end. 


Fig.  3.  CM-CFAR  block  diagram 


Fig.  4.  Complex  S/T  CFAR 


II.  Spatial/Temporal  CFAR 

The  proposed  method  was  shown  in  Fig.  4.  The  mathematical  analysis  is  given  in  below. 

A.  Analysis 

Theoretical  analysis  of  the  proposed  new  method,  when  CA-CFAR  is  used  in  the  spatial  averaging 
branch,  is  presented  here.  Assuming  a  square-law  detector,  Gaussian  noise  and  Swerling  II  model  for 
target,  the  relations  between  false  alarm  probability  and  detection  probability  in  optimum  case  are  [5] 

nOO 

Pfa  =  P[y  >  yo\H0\  =  /  l/2n  *  exp{-x/2p)dx  =  exp(—y0/2p)  (1) 

Jyo 

/•OO 

Pd  =  P[y  >  yn\Hi\  =  /  l/2/z(l  +  S)  *  exp{-x/2p(l  +  S))dx  =  exp(-y0/2p,(l  +  S)),  (2) 

Jyo 

where  yo  is  threshold,  p  is  noise  power  for  each  cell,  and  S  is  target  signal  to  noise  power  ratio.  Hi 
and  H0  represent  hypotheses  relating  to  the  presence  or  absence  of  target  respectively.  In  the  CFAR 
processor,  these  probabilities  are  computed  on  a  window,  and  are  then  averaged  on  the  probability 
density  of  Z ,  the  random  variable  of  reference  samples  over  all  windows.  In  the  complex  method 
presented  here,  Z  is  composed  of  two  random  variables,  Z\  and  Z2,  which  are  the  samples  that  partic¬ 
ipate  in  spatial  and  temporal  CFAR  respectively.  Therefore  the  detection  and  false  alarm  probability 
functions  are 


Pfa  =  Ez1,zAp[y  >  (aczi  +  0kz2)\Ho]} 

=  Ezi,z-2  {faczi+pkz2  l/2pexp(-y/2p)dy} 

=  EZltZ2{exp(-(aczi  +  j3kz2)/2p). 

Assuming  a  Swerling  II  model  for  the  target,  all  samples  are  independent.  Therefore 

Pfa  =  EZl{exp~{^l}}  *  EZ2{exp~(02htiZ2)} 

=  MZl{f^}*M22{§}. 

Similarly,  the  detection  probability  is  given  by 

Pd  =  EZltZ2{P[y  >  (aczi  +  /3kz2)\H1}} 

=  {  2^(1 +S)  }  *  AIZ2  {  2M(q  +  5)  }  • 

Applying  the  equations  given  in  [5] ,  [6] ,  for  the  CA-CFAR  processor 

MZl(u)  =  (l  +  2p^)~N, 


(3) 


(4) 


(5) 


(6) 


Threshold 


where  N  is  the  number  of  cells  that  are  used  in  the  estimation.  For  CM-CFAR 


OO 

MZ2(u)  =  JJ[1  +  2/j,uw(l  -  w)1]-1.  (7) 

1=0 

The  equations  for  false  alarm  and  detection  probability  are  thus  obtained  as 


“  m, 

Pfa  =  {lit1  +  aCW (1  ~  W)Z]_1}  *  (1  + 

1=0  ' 


(8) 


Pd  =  {III1 


1=0 


1  +  S 


(i  -«ot  }*(i  + 


(3k 


N(1  +  S) 


-N 


(9) 


In  the  rest  of  this  discussion  samples  of  1,000  scans  will  be  considered.  To  compare  the  complex 
method  with  CA-CFAR  and  simple  CM-CFAR,  the  parameters  of  the  two  methods  were  determined 
separately.  For  this  purpose,  with  Pfa  =  le  —  8  and  N  =  15,  the  following  were  obtained  (c  is  obtained 
via  numeric  calculation) 


k  =  36.2182 

c  =  31.5945  (10) 

w  =  0.125. 

Note  that  the  effective  number  of  integrated  pulses  in  CM-CFAR  is  Nf  =  [7].  Therefore,  w  was 

chosen  in  such  a  way  that  both  branches  in  the  complex  CFAR  (i.e. ,  CM-CFAR  and  CA-CFAR)  have 
the  same  integration  window  length. 

In  the  next  step,  a  and  (3  were  chosen  as  a  =  0.4  and  (3  =  0.6,  and  false  alarm  and  detection 
probabilities  were  computed  by  Eqs.  (8)  and  (9).  In  this  way,  it  can  be  observed  that  the  false  alarm 
probability  decreases  dramatically  and  reaches  1.1708e-10,  a  reduction  by  a  factor  of  85.41.  Fig.  5 
shows  the  curve  of  false  alarm  probability  versus  a  for  the  above  mentioned  parameters.  As  observed, 
the  false  alarm  rate  has  its  minimum  value  in  a  =  0.52  where  it  reaches  1/106.6  its  primary  value. 
Detection  probability  curves  for  three  methods,  CA-CFAR,  CM-CFAR  and  complex  CFAR  are  shown 
in  Fig.  6.  It  is  observed  that  the  detection  probability  curve  for  complex  CFAR  is  located  between 
two  other  curves,  therefore  the  detection  loss  of  complex  CFAR  is  less  than  CM-CFAR. 
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Fig.  5.  Pfa  versus  a  for  complex  S/T  CFAR 


Fig.  6.  Detection  probability  for  three  methods  in 
homogenous  environment 


B.  Self-Masking 

In  general,  the  density  function  of  a  LTI  filter  cannot  be  determined.  But  in  the  clutter-map  filter 
shown  in  Fig.  3,  the  central  limit  theorem  can  apply  when  w  goes  to  0,  because  in  this  case,  the 
integration  modifies  to  averaging  of  an  infinite  number  of  random  variables.  Therefore,  the  output 
density  function,  regardless  of  the  input  density  function,  has  a  Gaussian  distribution.  The  mean  and 
variance  of  the  filter  output,  as  shown  in  [7]  can  be  written  based  on  input  sequence  x(i)  as 


(11) 


n—  1 

Z(n  -  1)  =  w  ^(1  -  wy-^xti). 

i= 0 

Therefore  the  mean  and  variance  of  Z(n)  will  be 


Hz(n-1)  =w\J£(j,x  b=(  l-w)  (12) 

1  —  h2n 

al(n-l)=w2Y-wa2x.  (13) 

In  the  above  equations,  ax  and  \xx  are  input  sequence  statistics  of  filter.  When  n  — >  oo,  the  steady 
state  values  of  mean  and  variance  will  be 

M>z{n  -  1)  — »  Hz  =  Hx  (14) 

az(n  —  1)  *  &z  =  ax/Nf  Nf=2-^.  (15) 


To  illustrate  the  system  performance  in  presence  of  a  slow  target,  it  is  assumed  that  the  target  enters 
the  map  cell  in  the  nth  scan  and  stays  there  for  several  scans.  Therefore  the  filter  output  depends  on 
two  factors,  Background  estimation  in  steady  state,  and  the  number  L  of  scans  in  which  target  stays 
in  the  same  cell.  As  mentioned  in  [7],  in  this  case,  Eq.  11  modifies  to 

n— 1  n+L— 1 

Z(n  +  L  -  1)  =  w  6n+i-1-<n(i)  +  w  &n+L_1_ia;(i),  (16) 

i= 0  i=n 

where  n(i)  is  interference  signal  to  the  filter  input  before  that  target  enters  into  the  cell.  The  system 
has  been  assumed  to  be  in  steady  state  condition  before  the  target  enters  the  cell.  Thus,  the  statistics 
of  the  filter  output  can  be  written  as 

Hz(n  +  L  -  1)  ^  nn  *  bL  +  /ix  *  (1  -  bL)  (17) 

az(n  +  T  —  1)  —  jfjb2L  +  ^(1  —  b2L).  (18) 

In  the  rest  of  discussion,  the  mean  of  interference  and  target  will  be  considered  as  zero.  In  [7],  to  study 
system  behavior  in  the  presence  of  a  self-masking  target,  it  was  assumed  that  the  density  function  of 
the  threshold  is  still  Gaussian.  Therefore 


T(n  +  L  —  1)  oc  a),  (19) 

with  this  assumption,  the  variance  of  the  interference  after  the  target  entered  into  map  cell,  in  the 
Lth  scan  in  the  input,  can  be  written  as  [7] 

a2n,(n  +  L  —  1)  =  £ b2L  +  |l(l  -  b2L).  (20) 

Signal-to-noise  ratio  for  the  scans  in  which  the  target  has  been  in  map  cell  is 


S'  = 


(21) 


Replacing  this  value  into  Eq.  (9),  for  the  scans  in  which  target  is  present  in  the  map  cell,  gives  the 
system  performance  for  a  slow  target 


L—l 


^  =  { n  i1  +  rfwd  -  -)']-1  lb1  +  rbia  -  -T1}  a  + 

1=0  1  l=L 


(3k 


N(1  +  S) 


-N 


(22) 


where 


(23) 
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Fig.  7.  Detection  performance  for  a  self-masking  tar¬ 
get  for  SNR=40dB,  a  =  0.4 


Fig.  8.  Detection  performance  for  a  self-masking  tar¬ 
get  for  different  value  of  a ,  SNR=40dB 


Setting  (3  =  0,  the  system  performance  while  using  just  CM-CFAR  is  obtained.  Detection  probability 
curve  for  SNR=40dB  is  shown  in  Fig.  7.  As  is  evident  from  this  figure,  the  system  performance  for  a 
slow  target  in  complex  CFAR  is  better  than  CM-CFAR. 

Fig.  8  demonstrates  the  self-masking  effect  for  different  values  of  a  (0  <  a  <  1  with  step  0.1).  As 
expected,  as  CM-CFAR  use  increases  in  determining  the  threshold,  so  will  the  intensity  of  the  self- 
masking  effect.  Considering  Fig.  5  and  8,  the  parameter  a  can  be  determined.  Since  for  a  given  Pfa 
two  different  values  for  a  will  be  obtained,  to  decrease  self-masking  effect,  smaller  value  has  to  be 
chosen. 


III.  Experimental  Results 

Practical  observation  of  the  proposed  method  showed  that  by  using  just  CM-CFAR  (a  =  1, (3  =  0), 
slow  targets  like  helicopters  are  not  observed.  But  in  the  combined  method,  this  drawback  has  been 
diminished;  furthermore  ,  sub-clutter  visibility  has  been  increased.  In  this  observation  the  detection 
performance  of  the  proposed  method  was  not  less  than  CA-CFAR. 

IV.  Conclusions 

In  this  paper,  a  combined  method  of  spatial  and  temporal  CFAR  was  introduced  to  benefit  from 
their  advantages  at  the  same  time.  By  using  this  method,  three  advantages  are  clearly  evident: 
self-masking  problem  is  decreased;  sub-clutter  visibility  of  CM-CFAR  remains;  and  due  to  the  sliding 
window  in  the  spatial  branch,  periodic  response  is  eliminated.  While  using  spatial  CFAR  by  itself  sub¬ 
clutter  visibility  is  extremely  reduced;  however,  combining  two  methods  gains  both  of  their  advantages 
simultaneously. 
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